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The excitation curve of e + e~ annihilation into four charged pions in the p(770) region is calculated 
using three existing models with p mesons and pions in intermediate states supplemented by Feyn- 
man diagrams with the ai(1260)7r intermediate states. A two-term phenomenological Lagrangian 
of the aip-K interaction is used. The mixing angle is determined by fitting the e + e~ —* 7r + 7r~7r + 7r~ 
cross section data of the Novosibirsk CMD-2 collaboration and also its combination with the low- 
energy part of the BaBar collaboration data. It is shown that the inclusion of the am interme- 
diate states succeeds in obtaining a good agreement with the data on both cross section and the 
p° —¥ tt + tt~tt + tt~ decay width. When moving to energies above 1 GeV, the p(1450) and p(1700) 
resonances are taken into account to get excellent agreement with the BaBar data over the full 
energy range up to 4.5 GeV. 

PACS numbers: 13.30.Eg, 13.66.Bc, 12.39.Fe, 13.25.Jx 



I. INTRODUCTION 

Task of describing the excitation curve of the e + e~ — > ir + ir~ir + ir~ reaction at low energies is closely related 
to the investigation of the energy dependent decay width of the four pion decay of p(770). The validity of the 
factorization of the cross section into the p(770) production and decay parts is generally assumed on the basis 
of the vector meson dominance (VMD) hypothesis [1]. The assumption of factorization allows experimentalists 
to determine a particular partial decay width at the nominal p(770) mass on the basis of measurement the e + e~ 
annihilation cross section into the corresponding final state at the corresponding energy. The current experimental 
value r(p° — > 7r + 7r~7r + 7r~) = (2.8 ± 1.4 ± 0.5) keV was obtained in that way [2]. In this work, we will also assume 
a one-to-one correspondence between the cross section and the energy dependent decay width at the same energy, 
ignoring the complications which may appear if some conditions are not met [3] . Needless to say that the evaluation 
of the decay width is less demanding technically and computationally than that of the cross section (five-dimensional 
quadrature instead of eight-dimensional one in the case of four-body final states). 

The cross section of the e + e~ annihilation into the 27T+27T - and 7r + 7r~27r° final states was considered by Decker, 
Heiliger, Jonsson, and Finkemeier [4] in conjunction with the CVC-related decays of the r lepton. The intermediate 
states of their model contained p(770), ai(1260), and a scalar-isoscalar two-pion resonance. In the two-charged-two- 
neutral case also the w(782) was included. The aip7r vertex factors were adopted from the study of the three-pion 
decay of the r lepton by Isgur, Morningstar, and Reader [5]. Czyz and Kiihn [6] constructed a Monte Carlo generator 
of the reaction e + e~ — > 7 + 4-7T. The hadronic matrix elements were used in the form suggested in [4], corrected 
only for some minor deficiencies. To check the soundness of their approach, Czyz and Kiihn calculated also the 
excitation curves of the nonradiative reactions e + e~ — * 7r + 7r~7r + 7r~ and e + e~ — ► 7r + 7r~27r° in qualitative agreement 
with available data. Ecker and Unterdorfer [7, 8] performed the first calculation of the processes e + e~ — > 47r and 
t — * v t \k with the correct structure to 0(p 4 ) in the low energy expansion of the Standard Model extrapolated to 
the resonance region. To get a good description of the e + e~ — > 7r + 7r~7r + 7r~ cross section up to 1 GeV, they had to 
include as additional contribution the a\ exchange. They circumvented the aip7r Lagrangian ambiguity by choosing 
a special relation among the individual coupling constants. 

The four-pion decays of the p(770) are generally considered a convenient test ground of the low-energy effective the- 
ories of the interactions of p mesons and pions. In the past, several papers appeared that calculated the corresponding 
partial decay widths [9-13]. Moreover, Achasov and Kozhevnikov [12, 13] argued that the four-pion decay widths 
of p(770) are not experimentally well defined because they require the averaging over a mass interval in which they 
rise rapidly. They therefore calculated, in addition to the decay widths r(p° — > 7r + 7r~7r + 7r~), r(p° — > 7r + 7r~7r 7r°), 
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T(p + — > 7r + 37r ), and — > 27r + 7r~7r°), the e + e~ — > 7r + 7r~7r + 7r _ reaction cross section as a function of incident 
energy (excitation curve) and compared it to the CMD-2 data [2] from Novosibirsk. 

With respect to the role of the axial- vector meson ai(1260) in the four-pion decays of p°, the situation is somewhat 
controversial. On one side, the intermediate states containing the ai(1260) were either ignored [9-13] or shown [11] 
to have little influence on the four-pion decay widths of p(770) 1 . On the other side, the analysis of the differential 
distributions of charged pions coming from the e + e~ annihilation in the energy range 1.05-1.38 GeV demonstrated 
the dominance of the ai(1260)7r intermediate states [14]. Given the large width of ai(1260) (T ai =250 to 600 MeV 
[15]) it would be surprising if the role of the ai-meson diminished so fast outside the above range, which is not too 
far from the p(770) mass. Moreover, the ai(1260) meson was shown to be important in the four-pion decays of the 
r-lepton [4, 16, 17], which are in a sense isospin counterparts of the four-pion final states in the e + e~ annihilation. 
Recently, the paper of Achasov and Kozhevnikov has appeared [18] that included the intermediate states with the a\ 
meson using the generalized hidden local symmetry model [19]. This increased the p° — > 2ir + 2ir~ decay width from 
0.94 keV to 1.59 keV assuming the nominal mass of the a± resonance, see Table I in [18]. Unfortunately, the authors 
of [18] do not provide the comparison with the e + e~ — > 2ir + 2ir~ cross section, which has a greater discriminatory 
value than the decay width alone [12, 13]. 

The present work was triggered by our interest in the electromagnetic probes in relativistic heavy- ion collisions. 
The production of prompt dileptons and photons is for a long time considered a powerful tool for investigating 
the properties of dense systems created in the hadronic and nuclear collisions [20, 21]. The interest of the heavy 
ion community in dilepton production has recently been boosted by very precise dimuon data by the CERN/NA60 
collaboration [22]. The prompt dileptons and photons can, in principle, originate from two sources: (i) quark-gluon 
plasma, (ii) hadron gas. The theoretical calculations of the dilepton and photon yield from the latter are hampered 
by not uniquely known Lagrangian of the aipn interaction [23, 24]. We suggest one way how to relieve this problem 
and make the predictions of the dilepton and photon production from hadron gas more reliable. The results of the 
present work have already been utilized [25] in the evaluation of the dimuon production rate in In-In collisions at 158 
A GeV. 

A way of narrowing the uncertainty interval in dilepton production from the hadron gas was advocated a long 
time ago by Li and Gale [26]. They checked the feasibility of the dilepton rate calculation in, e.g., ojtt collisions, by 
comparing the inverse process (e + e~ annihilation into the uitt° final state) with the available experimental data. Li 
and Gale also considered the e + e~ annihilation into four pions in the narrow a\ width approximation, i.e. as process 
e+e~ — > a\w. In the present work we handle it as a genuine four-fmal-pion process, considering not only the a\TT 
intermediate states but also other intermediate states following from various models based on the chiral perturbation 
theory. 

In this paper we investigate the role of the ai(1260) resonance in the e + e~ annihilation into four charged pions and 
in the four-charged-pion decays of p(770) in more detail. We supplement three existing models, which consider only 
p and 7r in intermediate states (diagrams (a) and (b) in Fig. 1), with the ai contribution (diagrams (d) in Fig. 1). 
Those three models are (1) the model of Eidelman, Silagadze, and Kuraev [10], (2) one of the models considered by 
Plant and Birse [11], and (3) the model of Achasov and Kozhevnikov [12, 13]. We will consider only the final state 
with all charged pions, for which the experimental data are best. The difference between the approach of Achasov 
and Kozhevnikov [18] and ours lies mainly in the aipn Lagrangian. In [18], the generalized hidden local symmetry 
Lagrangian was chosen, whereas we choose a more phenomenological two-term Lagrangian, which is often used in a 
different branch of the particle physics. Its individual terms and their specific combinations appeared in many papers 
computing the dilepton and photon production rate from a thermalized meson gas. In this paper we consider the 
mixing angle between the two terms as a free parameter and will determine its value by fitting the excitation curve 
of the e+e~ — > 7r + 7r~7r + 7r~ reaction. 

In order to evaluate the amplitude induced by eight diagrams Fig. 1(d) we have to choose a Lagrangian of the aipir 
interaction. For this choice, there are basically two approaches in the literature. One explores well defined theoretical 
concepts to build dynamical models, the free parameters of which are then fixed by comparison with observed masses 
and decay widths. See, e.g., Refs. [23, 27-34]. 

In other, more phenomenological, approaches the authors simply chose for the ctipn Lagrangian various expressions 
built from the field operators and compatible with the fundamental conservation laws. Such Lagrangians, after 
fixing their coupling constants, were then used to calculate various observable quantities, see, e.g., [35-38]. In some 
approaches [5, 39], directly the vertex factors were written without showing the corresponding Lagrangian. From the 
fact that different authors pick different Lagrangians one might get the impression that the choice of Lagrangian is 
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FIG. 1: Selected Feynman diagrams describing decay p° — > ty + it~ ty + tt~ 

not very important and that various Lagrangians lead to identical, or at least similar, results. This is not true, and 
the observable quantities may be very sensitive to the choice of the aipir Lagrangian, as was demonstrated, e.g., by 
Song [23]. The discussion of the a\ phenomenology with emphasis on the photon and dilepton production from a hot 
meson gas can be found in Ref. [24]. 

With so many aipir effective Lagrangians, it is interesting to learn which Lagrangian is preferred experimentally. 
This would require constructing a general Lagrangian with a set of free parameters and fixing them by comparing all 
possible observables with the existing data. Of course, such a program is very ambitious. In this paper we are going 
to do something much simpler. Below, we choose a two-component Lagrangian and determine its two free parameters 
by requiring that the decay width of ai(1260) be reproduced and the best possible fit obtained for the excitation 
curve of the e + e~ — > ir + ir^ir + ir~ reaction. Even this restricted program cannot be accomplished completely. Firstly, 
the width of ai(1260) is not known reliably. We will consider three values from the range given in the Particle Data 
Group tables [15], namely 250, 400, and 600 MeV. Secondly, the result of the fit will depend also on the basic p and 7r 
intermediate state model to which we add the aiir contribution. Nevertheless, we will show that the inclusion of the 
ai7r intermediate states is necessary for obtaining good agreement with the e+e~ — > ir + ir~ir + ir~ excitation curve. 

Following the suggestion of T. Barnes [40] we also include the ratio of the Z?-wave and S'-wave amplitudes of the 
a\ — > pir decay as a fitted quantity. The importance of the D/S ratio for selecting among the aipir Lagrangians was 
stressed in a different context in [24]. We use the experimental value D/S = —0.14 ±0.11 [41], which was obtained by 
genuine partial wave analysis of the reaction n~p — ► 7r + 7r~7r~p. We consider the other values which exist in literature 
strongly model biased. They were obtained by fitting the three-pion mass spectrum in the decay t — ► ^ T 37r using the 
model of Isgur, Morningstar, and Reader [5] and then calculating the D/S ratio from the optimal parameters of the 
model. 



II. ORIGINAL MODELS, MODIFICATIONS, AND ADDITIONS 

As we already stated, we will complement three existing models of the four-pion decays, which consider only 
intermediate states with p mesons and pions, with the intermediate states containing the axial- vector resonance 
ai(1260). We must admit that our choice of models is rather arbitrary. Moreover, we took the models as they 
appeared in the literature and did not check their compatibility with chiral symmetry and with the constraints on 
coupling of resonances to the pions [42, 43]. 

Those three models are characterized below. 
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A. Model of Eidelman, Silagadze, and Kuraev (ESK) 

This model [10] is based on the effective chiral Lagrangian by Brihaye, Pak, and Rossi [44], which was investigated 
also in [45]. It follows from that Lagrangian that all (a) and (b) diagrams depicted on Fig. 1 contribute to the 
p° — ► 7r + 7T~7r + 7r~ decay rate. Their amplitudes (in the notation slightly different from ours) are shown in the paper. 
Our usage of this model will differ from the original paper in three respects: (1) We add a\ diagrams Fig. 1(d). (2) 
We use a different value of the parameter ah, defined in [10]. Instead of 0.55 we set = 0.5, which follows from 
the KSRF relation [46], to be in conformity with other two models. (3) We replace the scalar part of the p-meson 
propagator with fixed mass and fixed width by the prescription 

Pp{s) = ' s -M^ s ) + imp T p ( s y {2A) 

which uses the running mass squared M 2 (s) and the energy dependent total width T p (s) from Ref. [47]. 

The last point deserves more comments. The denominator of our propagator (2.1) is an analytic function in the 
s-plane with a cut running from Arn^ to infinity, as required by general principles. This property differs (2.1) from 
most of the formulas used in the literature. The real function M 2 {s) is calculated from T p (s) using a once subtracted 
dispersion relation, which guarantees that the condition M 2 (m 2 ) = m 2 is satisfied. Further condition 



dM 2 (s) 



ds 



(2.2) 



is not fulfilled automatically and serves as a test that all important contributions to the total p-meson width T p (s) 
have properly been taken into account. See [47] for details. If we replace the m p accompanying T p (s) in Eq. (2.1) by 
■y/s, as it is done in some existing formulas, the condition (2.2) cannot be satisfied for any reasonable choice of T p (s). 

The running mass approach [47] takes into account, in addition to the basic two-pion decay channel, several channels 
(w7T°, K + K~, K°K°, and i]ir + Tr~) which open as the p resonance goes above its nominal mass. It also considers 
structure effects described by the strong form factors. In these two respects it differs from other approaches that 
appeared in the literature [48-50]. Gounaris and Sakurai [48] considered only the two-pion contribution to the total 
width of the p a resonance and ignored structure effects. Vaughn and Wali [49] took into account the strong form 
factor, but again ignored higher decay channels. Melikhov, Nachtmann, Nikonov, and Paulus [50] included the K + K~ 
and K°K° channels, but did not consider the strong form factors. 

We will use the p propagator (2.1) not only in conjunction with the ESK model, but also with the other ones. This 
is the main reason why our results calculated within the original models (i.e., without the aiir intermediate states) 
and presented below differ slightly from the results quoted in the original papers. 



B. One of the models of Plant and Birse (PB/HG) 

Plant and Birse [11] investigated several models of the four-pion decays of p°(770). One of them (labeled HG) is a 
corrected version of the model by Bramon, Grau, and Pancheri [9], which was based on the hidden gauge theory of 
Bando et al. [51]. The 2ir2p contact terms, see diagram (bl) in Fig. 1, are missing in this approach. The amplitude 
of the (al) diagram is different from that in work by Eidelman, Silagadze, and Kuraev [10] by a factor (—1/2). The 
amplitudes of (a2) and (b2) diagrams are equal to their ESK counterparts. 



C. Model of Achasov and Kozhevnikov (AK) 

Achasov and Kozhevnikov [12, 13] studied the four-pion decays of p(770), five-pion decays of w(782), and the 
processes related to them. Namely, the e + e~ annihilation into the four- and five-pion final states and the four-pion 
decays of the r lepton. They used the Weinberg Lagrangian [52] obtained upon the nonlinear realization of chiral 
symmetry. From their rather extensive work we adopt their prescriptions for the amplitudes (al), (a2), and (b2) 
of the p° — > 7r + 7r~7r + 7r~ decay. The contact amplitude (bl) is again vanishing. Achasov and Kozhevnikov used a 
fixed-mass, variable-width formula for the p-meson propagator. 
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D. a\p-K Lagrangian and the amplitude of the diagrams containing the a\ meson 

We choose the following interaction among the a\, p, and tt fields 

C=^^(Cicos6 + C 2 sm9), (2.3) 

where .g aip7r and 9 are yet undetermined parameters, 

A = A" • (V„„ x 9» , (2.4) 
C 2 = V M „ • (d»A» x 0) , (2.5) 

and V M „ = <9 M V„ — d v V p . The isovector composed of the p-meson field operators is denoted by V M , similar objects 
for 7T and ai are and A**, respectively. We write 

01 = ^|(0c + 0c), 

02 = ^ (0c - C ) , 

03 = 0n, 

and assume that C contains the annihilation operators of the positive pion and creation operators of the negative 
pion. <p n is the operator of neutral pion field. 

A specific combination of terms (2.4) and (2.5) appeared in the pioneering work by Wess and Zumino [27]. Term 
(2.4) alone was used by Xiong, Shuryak, and Brown [39] in their study of the photon production from meson gas. 
Janssen, Holinde, and Speth [37] picked the term (2.5) when they evaluated the amplitude of the Trp scattering. 
Another combination of (2.4) and (2.5) appeared in the calculation of dilcpton production from meson gas by Song, 
Ko, and Gale [53]. 

Lagrangian (2.3) leads to the following factor for the vertex in which an incoming (index a), an outgoing p° 
(index p) , and an outgoing 7r + meet 

V a »(p ai ,p p ,p„) = i^{ cos e[p a p p^-(p n p p )g^] 



— sin i 



The a l p°n vertex acquires an extra minus sign. The evaluation of the decay rate of af — > p°ir + using vertex (2.6) 
is straightforward. 



2 



°i w+ 1927rm3 



ai 



where 



and 



x R{m 2 ai ,m 2 p ,ml+) , (2.7) 
\(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz 



R(x,y,z) = \{x-y- zf + — (ar-j/ + zf 



cos 2 I 



2 [(x - z) 2 + y(x + z - 2y)] cos6>sin6> 
+ [(x + y- z) 2 + 2xy] sin 2 9 . 

If we assume the charge independent aip7r coupling constant and masses of p and ai, the width of the decay — > p + ir a 
is obtained from (2.7) by changing just the pion mass. Formula 

r a+ = r a+^p°7r+ + r a+^p+ TrO ( 2 - 8 ) 

enables us to find the coupling constant g aiplT for given T a + and sin0. Because each of the diagrams Fig. 1(d) contains 
two aipK vertices, the overall sign of the Lagrangian (2.3) is not important and we can assume a non-negative cos#. 
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A question arises whether the narrow p-width approximation used above is accurate enough for the purpose of 
determination the g aip7T coupling constant. Achasov and Kozhevnikov [18] showed that the a\ decay width calculated 
as r(ai — ► pn) came out larger than that calculated as T(ai — ► 37r) using the same coupling constant g aiP -n- We have 
examined this issue, too, using our two-component Lagrangian (2.3) and got that the T{a\ — > 37r)/r(oi — > pir) ratio 
is smaller than unity (like in [18]) for sin# < 0.2 and sin^ > 0.65. In the remaining interval, which contains all the 
values of smO that will be met in our calculations, this ratio is greater than one. Moreover, we have found that if one 
takes into account the strong form factors in the aipn and pirn vertices, the results of both approaches become almost 
identical. This can be explained as follows. According to Kokoski and Isgur [54] formula, which will be shown later, 
the strong form factor in a particular decay vertex is a decreasing function of the three-momentum of an outgoing 
particle in the rest frame of the parent particle. In the <Zi — > 37r decay the intermediate mass of p is mostly smaller 
than its nominal mass, what means higher momenta of particles emerging from the aipn vertex. In addition, the 
a\ — > 3tt width is further reduced by the form factor in the pirir vertex. 

Let us now turn to the amplitude of the a\ diagrams Fig. 1(d) for the decay p°(p) — > n~ (pi)n + {p2)^ + (P3)^~ {pa) ■ 
We first introduce the notation 

Qi = P-Pi, 

T*ij Pi ~\~ Pj , 

S ij = r iji (2-9) 



and then write the amplitude in the form 

( (A) 

where is the polarization vector of the decaying p° and 



Jd,» = -(l-i , 12^4)(l + i , u)(l + P23)V r aM (-«4,-p,P4) 

x P a f(94)V(94,r 12 ,p 3 ) 3p (p 2 - Pl yp P (s 12 ). 

Here, Pij denotes the operator that interchange four-momenta pi and pj. The axial- vector meson propagator 

Pafil) = l-o ~ ^- (2-10) 

is chosen in a simple fixed-mass, fixed-width form. Here, we are going to consider the four-pion system with invariant 
energies less than 1 GeV. The invariant mass of the three-pion system, which is equal to \fq* \ is thus limited by 
0.86 GeV. Achasov and Kozhevnikov [18] showed that the ai decay width is negligible in that energy range. Referring 
to their finding we set T ai = in Eq. (2.10). The scalar part of the p propagator is again used in the form (2.1). 

E. Technicalities 

The complete amplitude of the p a — > 7r + 7r~7r + 7r~ decay is 

M^=e^, (2.11) 
where is the polarization vector of the decaying p° and 

Four-vectors J a ^ and J^^ describe the contributions from (a) and (b) diagrams in a particular model and Jd, P is the 
contribution from (d) diagrams. The sum over the p-meson polarizations of the amplitude (2.11) squared is given by 



2 / p VpV 



^ + ^)«- (2-12) 



This formula is more complicated than that used in [10], because the four- vectors J a41 and J 6jA1 of PB/HG and AK 
models do not satisfy the transversality condition J^p^ 1 = 0. We used the algebraic manipulation program REDUCE 
[55] to express the sum (2.12) in terms of six invariants s^-, i < j, j — 2,3,4 defined in (2.9). Of course, only five of 



7 



them are independent and we used the identity Y^ii<j s ij = m p + 877?,^ for the checks in the process of evaluation of 
the decay width. When calculating the excitation curve, m 2 p is replaced by s, the square of the incident energy. 

When calculating the decay width of an unpolarized parent particle, we may take advantage of the spherical 
symmetry of the problem and choose the following kinematic configuration: (1) The parent particle a is at rest. (2) 
The summed momentum pi2 of particles 1 and 2 points in the direction of the z axis. (3) The individual momenta 
Pi and P2 lie in the xz plane. Then the following formula, written in a general case with arbitrary masses and spins, 
is valid 



N 



dm\2 Pi / dm 3 4 



16(2"7r) 6 m^ Jm 1 +m 2 Jm.3+m 4 

dcosdl / dcostfg / dip 3 \M\ 2 . (2.13) 
-1 J-l Jo 

The last quantity is the amplitude squared, averaged over the initial spin states, and summed over the final spin 
states. The factor N takes into account the identity of the final particles and equals 1/4 in our case. The asterisk 
denotes the momentum in the corresponding rest frame (1-2 or 3-4), p — |p|, and rriij — ^/sjj = yj (pi + Pj) 2 - 

For evaluation of the integrals in (2.13) we used a sequence of the five one-dimensional Gauss-Legendre quadratures 
of the sixteenth order. We prefer this method to Monte-Carlo integration because we use the result of the integration 
in a minimization procedure and therefore we require that the same value of the optimized variable (Lagrangian mixing 
angle) yield always the same value of the minimized function, which would not be satisfied with the Monte-Carlo 
integration. Nevertheless, we checked our computer code by evaluating the decay width for a particular value of the 
mixing angle using a completely independent code based on the Monte-Carlo method. 

To convert the calculated decay width into the cross section, we start with the formula 



o-47r(s) = 7rr 4lr (s). 

1 TT + TT- (S) 



Using 



3/2 

2 



where F 7T (s) is the contribution of the p resonance to the pion form factor, and 

3/2 



g 2 p W / 4m»Y 



with W = y/s, we arrive at 



We further use the VMD expression for the dielectron decay width of p° 



4nm ( a x 



and get 



- ^^tr WAs)\ 2 T^(s). (2.16) 



TTlpW 3 



If we set, following Achasov and Kozhevnikov [13], 



2 

m„ 



FAs) - j^ry (2.17) 



where the inverse p-meson propagator D p (s) is defined in Eq. (2.2) of [13], we reproduce their Eq. (3.1). 
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TABLE I: x 2 /NDF of the fits to the CMD-2 cross section data (11 data points) 



r 

1 ai 


L/O-LV 


PB /HG 


AK 


uiiiy LL± 


(MeV) 


[10] 


[11] 


[12, 13] 




250 


1.60 


1.34 


1.28 


1.68 


400 


1.53 


1.37 


1.30 


1.82 


600 


1.61 


1.41 


1.31 


1.94 


Only p, TV 


17.6 


15.0 


14.8 


/ 



In our opinion, Eq. (2.16) overestimates the cross section if the experimental value of r e + e - is used. The reason is 
that the dielectron decay width calculated from (2.15) is smaller than the experimental value. We will therefore stick 
with formula (2.14). 

We utilize the scalar part of the p-meson propagator (2.1) to write our Ansatz for the p- meson contribution to the 
pion form factor 

M 2 (0) 

FAS) = M?(s)-s-im p T p (sy (2 - 18) 

As shown in [47], this formula gives the correct value of the mean square radius of the pion. The form factor (2.17) 
fails in this test. 

For the p coupling constant we use the same value as in [11-13], namely g p — 5.89. This value is compatible 
with what follows from the KSRF relation [46] (5.900 ± 0.011). Both values are little lower than g p = 6.002 ± 0.015 
calculated from the p-meson width. 



III. LOW-ENERGY RESULTS (W < 1 GeV) 

We deal with the excitation curves of the reaction e+e~ — ► 7r + 7r~7r + 7r~ calculated in three different models (ESK, 
PB/HG, AK) supplemented with the a\ diagrams Fig. 1(d). We first fit them to the CMD-2 data [2] by varying the 
sine of the mixing angle 8, defined in (2.3), for the three fixed values of the width of the ai(1260) meson. We did 
not consider the first two points in the CMD-2 data, because they give only upper bounds of the cross section. The 
ratios of the usually defined \ 2 to the number of degrees of freedom (NDF), which characterize the quality of the 
fit, are shown in Table I. The last row in Table I shows the values of x 2 /NDF that indicate how well (or badly) the 
original models without a\ agree with the data. No free parameter is involved in the latter case. Table I shows that 
the ratio x 2 /NDF is always greater than one. In what follows, we shall therefore multiply the statistical errors of the 
quantities obtained in the process of minimization by the square root of that ratio. 

The inspection of Table I shows that the inclusion of the ai contribution greatly improves the agreement with the 
data. The interference between the original diagrams and the new ones is important, the results of the combined 
model are better than those of the a\ diagrams alone. The best results (lowest x 2 ) are obtained with the AK model 
supplemented with the aiu intermediate states. 

To investigate the sensitivity of our results to the input data, we combine the CMD-2 data [2] and the low-energy 
(s < 1 GeV 2 ) part of the BaBar data [56] (BaBar-LE in what follows) into a new set and repeat the calculations. 
The results are shown in Table II. Their comparison with the results obtained from the CMD-2 data alone shows 
two important differences: (i) The agreement of all models with data, characterized by x 2 /NDF, is now better. It 
indicates that the CMD-2 and BaBar-LE data are compatible, so the increased number of data points does not bring 
proportional increase of x 2 - (ii) Whereas merging of the a\ diagrams with the PB/HG or AK model improves the fit, 
adding the ESK model diagrams to the pure a\ contribution leads to the opposite effect. 

Next, we add the D/S ratio [41] to the set of fitted experimental values and repeat the calculations. The results are 
shown in Tabic III. The salient feature of those results is a clear preference of the highest assumed value (600 MeV) 
of the total a\ width. 

In Figs. 2,3, and 4 we show the comparison of the combined set of data with the excitation curves calculated in 
all three models supplemented with the a\ diagrams Fig. 1(d). The same comparison for the a\ diagrams alone is 
depicted in Fig. 5. 

In Table IV we can see the values of sin 6* together with their errors (defined in the usual way [57]) obtained from 
the fit to the CMD-2 & BaBar-LE cross section data and to the D/S ratio. As we mentioned above, three different 
values of the ai(1260) width are assumed. Table V compares the experimental value of the p° — ► 7r + 7r~7r + 7r~ decay 
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TABLE II: x 2 /NDF of the fits to the combined CMD-2 & BaBar-LE cross section data (27 data points) 



r 01 






AK 


only oi 


(MeV) 


[10] 


[ii] 


[12, 13] 




250 


1.42 


1.20 


1.19 


1.32 


400 


1.48 


1.21 


1.18 


1.39 


600 


1.55 


1.22 


1.19 


1.44 


Only p, TV 


9.4 


10.4 


10.3 


/ 


TABLE III: xV NDF of the fits to the CMD-2 & BaBar-LE cross 


section data and to the D/S ratio (28 data points) 




ESK 


PB/HG 


AK 


only ai 


(MeV) 


[10] 


[11] 


[12, 13] 




250 


3.66 


1.95 


1.99 


1.96 


400 


1.98 


1.34 


1.33 


1.48 


600 


1.65 


1.20 


1.18 


1.41 



width [2] with the results obtained from various models under the same conditions. The results for both quantities 
obtained with other data sets (CMD-2 only, CMD-2 & BaBar-LE) are very similar. 

IV. HIGH-ENERGY RESULTS (W UP TO 4.5 GeV) 

When we want to get a good description of the data on the e + e~ annihilation to four charged pions at energies 
above 1 GeV, we should consider also the contribution from diagrams where higher p resonances couple to the virtual 
photon and then convert into four pions. We include two resonances: p' = p(1450) and p" = p(1700). We assume 
that the decay of those resonances into four pions is governed by the same Feynman diagrams as that of p(770), with 
all coupling constants scaled by the same factor (different for p' and p"). This assumption implies that the four-pion 



1 . 1 M i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i ij 




W (GeV) 

FIG. 2: Excitation curves calculated in the original (without ai meson) and expanded ESK model compared to the CMD-2 
and BaBar-LE data. The D/S ratio was also used in fit. 
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0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 
W (GeV) 



FIG. 3: Excitation curves calculated in the original (without ai meson; dash-dotted curve close to the abscissa) and expanded 
PB/HG model compared to the CMD-2 and BaBar-LE data. The D/S ratio was also used in fit. 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 
W (GeV) 



FIG. 4: Excitation curves calculated in the original (without a\ meson; dash-dotted curve close to the abscissa) and expanded 
AK model compared to the CMD-2 and BaBar-LE data. The D/S ratio was also used in fit. 



TABLE IV: Values of sin 9 from the fit to the CMD-2 & BaBar-LE data and to the D/S ratio. 





ESK 


PB/HG 


AK 


only ai 


(MeV) 


[10] 


[11] 


[12, 13] 




250 


0.4092(33) 


0.4278(32) 


0.4267(32) 


0.4312(35) 


400 


0.4352(24) 


0.4624(34) 


0.4608(32) 


0.4679(39) 


600 


0.4659(27) 


0.5046(44) 


0.5022(41) 


0.5132(55) 



11 



1 . 1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



1 r e + e -> tt + tt 7t + 7t 




W (GeV) 

FIG. 5: Excitation curves calculated from the aiir diagrams only, compared to the CMD-2 and BaBar-LE data. The D/S ratio 
was also used in fit. 



TABLE V: Decay width r(p° — > tt + ty ti + ty ) (keV) calculated in various models using sin# from the fits to the CMD-2 Sz 
BaBar-LE data and to the D/S ratio. Experimental value is (2.8 ± 1.4 ± 0.5) keV [2]. 



r B1 


ESK 


PB/HG 


AK 


only ai 


(MeV) 


[10] 


[11] 


[12, 13] 




250 


4.28(01) 


3.16(25) 


2.70(23) 


4.52(30) 


400 


2.81(01) 


3.55(28) 


3.03(26) 


5.08(32) 


600 


1.94(02) 


3.77(30) 


3.22(27) 


5.39(37) 


Only p, TV 


16.2 


0.59 


0.89 


/ 



decay widths of p' and p" have the same shape in W as that of p(770). They only differ from it by constant factors. 
The simplifying assumption we have made allows us to use the same cross section formula (2.14) as in the low energy 
case, with F 7T (s) replaced by 

F(a) = F„(a) + 6F p ,(s) + eF p „(s), (4.1) 
where F p (s) differs from (2.18) by including also T p o ( s ) into the total decay width of p(770), 

2 

m„, 

P,( s )= E (4.2) 

and similar expression holds also for F p n(s). Constants 5 and e not only include the coupling constants modification 
factors mentioned above, but also account for the fact that the couplings of p' and p" to photon differ from that of 
p(770), which is fixed by VMD. Unknown complex parameters 5 and e will be determined, together with the masses 
and widths of p' and p" and other two parameters mentioned later on, by fitting the experimental excitation function 2 . 



2 When we replaced (2.18) by (4.1) in the low-energy region and kept the form-factor parameters as determined in the high-energy fit, 
we got the results that differed slightly from that in Sec. III. If we varied also those parameters when fitting the low-energy data, they 
acquired unphysical values (masses of p' and p" around 1 GeV) . It may signify that some contribution important at low energies (scalar 
resonances) is still missing in our approach. 



12 



TABLE VI: Results of the fit to the BaBar cross section data and to the D/S ratio (145 data points) for Y ai = 600 MeV. 



Model 


ESK 


PB/HG 


AK 


only anr 




[10] 


[11] 


[12, 13] 




X 2 /NDF 


1.21 


1.12 


1.12 


1.12 


sin 8 


0.4474(22) 


0.4592(28) 


0.4588(27) 


0.4603(28) 


(3 (GeV) 


0.3505(89) 


0.3665(97) 


0.3657(97) 


0.3695(98) 


m p , (GeV) 


1.419(12) 


1.439(13) 


1.438(13) 


1.442(13) 


I> (GeV) 


0.564(20) 


0.568(21) 


0.568(21) 


0.566(21) 








n 1 1 44^1 


0. 1 liy ^OtD 1 


lm(S) 


-0.039(11) 


-0.019(12) 


-0.021(12) 


-0.015(13) 


rn p „ (GeV) 


1.903(21) 


1.923(24) 


1.922(24) 


1.926(24) 


F p „ (GeV) 


0.247(38) 


0.284(44) 


0.283(44) 


0.290(45) 


Re(e) 


-0.0016(11) 


-0.0002(17) 


-0.0003(17) 


0.0002(18) 


Im(e) 


-0.00373(94) 


-0.0054(12) 


-0.0054(12) 


-0.0056(13) 



Another effect that has to be taken into account when dealing with higher energies is connected with the structure of 
the strongly interacting particles. Our decay amplitudes have been derived under the assumption that the pions, p's, 
and ai's are elementary quanta of the corresponding quantum holds. But this assumption is justified only when their 
mutual interaction is soft. When the momenta of the mesons which enters a specific interaction vertex get higher, 
the contribution of that vertex to the amplitude becomes smaller than in the case of the point-like participants. 
This effect is usually described by strong form factors. Given the present status of the strong interaction theory, we 
have to turn to models. For example, in the chromoelectric flux-tube breaking model of Kokoski and Isgur [54], the 
vertex describing a two-body decay is modified by the factor exp{— p* 2 /(12/3 2 )}, where p* is the three-momentum 
magnitude of the decay products in the parent particle rest frame and (3 w 0.4 GeV. For decay of the p meson with a 
(non-nominal) mass W into two on-mass-shell pions, this form factor can be written as 

F *'( s ) = ex p{-^r}< ( 4 - 3 ) 

where s = W 2 and Sq is the threshold value of s (4m 2 in the two-pion decay). The complete amplitude of the four-pion 
decay of p° contains many vertices, some of them with more than three incoming/outgoing particles. Applying the 
Kokoski-Isgur factor to each of them would be cumbersome and would require additional assumptions in the case of 
more complicated vertices. We will therefore assign an "effective" strong form factor of the form (4.3) to the complete 
amplitude of the decay p° — > n + n~ tt + n~ , but with s = 16m 2 . When fitting the experimental excitation curve, (3 
will be considered as another parameter. With the masses and widths of p' and p" , complex parameters <5 and e, and 
with the sine of the mixing angle 8 there are ten real parameters to be determined by fitting the e + e~ — > 7r + 7r~7r + 7r~ 
cross section data of BaBar collaboration [56]. 

In the following, we assume the a\ decay width of 600 MeV, for which the results of all models at energies below 
1 GeV were best. The same value was used in the a\ propagator (2.10). 

The resulting optimized values of the parameters listed above and their MINUIT [57] errors are shown in Table VI for 
the three different models supplemented with the aiir intermediate states and for the latter alone. Mutual comparison 
of the % 2 /NDF ratios clearly shows that the presence of the aiir intermediate states is crucial for obtaining good 
agreement of the calculated excitation curve with data. They provide a good fit even if taken alone. Adding the 
diagrams with 7r's and p's in the intermediate states does not change the quality of the fit if their amplitudes are taken 
from the PB/HG and AK models. On the other hand, the inclusion of the amplitudes of the ESK model brings some 
deterioration of the fit. The values of parameter (3 do not differ very much from the value advocated in [54] for the 
three-line vertices, what indicates that the effective-strong-form-factor approach we have chosen (4.3) is reasonable. 
The values of the sine of the mixing angle 9 are somewhat lower than those at low energies. The central values of 
the masses and widths of p' and p" a little different from those listed in [15], but differences are acceptable keeping 
in mind relatively large errors. 

The graphical comparison of the excitation curve calculated from the model containing only the diagrams with the 
oi7r intermediate states, shown in Fig. 1(d), with experimental data [56] is presented in Fig. 6. The excitation curves 
of other models (ESK, PB/HG, AK) combined with the qitt contribution differ only slightly and are not shown. 



13 




1 1.5 2 2.5 3 3.5 4 4.5 
W (GeV) 



FIG. 6: Theoretical excitation curve compared with the BaBar data [56]. The D/S ratio was also used in fit. The result of the 
pure ai model is shown. The other models combined with ai7v intermediate states provide almost identical curves. 

V. CONCLUSIONS AND COMMENTS 

Our low-energy results show that the inclusion of the ciiir intermediate states is of vital importance for obtaining 
a good agreement with the experimental data on the cross section of the reaction e + e~ — > tt + tt~ tt + tt~ as a function 
of the incident energy (see Tables I, II, and III or Figs. 2, 3, and 4). The x 2 /NDF ratio gets much smaller if a 
particular model is supplemented with the diagrams containing the a\ resonance in the intermediate states. Viewing 
from another perspective, the pure ai model provides relatively good agreement with the cross section data, much 
better than each of the three models without the am intermediate states. See Tables mentioned above and Fig. 5. 
Adding the diagrams from the original models to the pure a\ model improves the fit in the case of the PB /HG and 
AK models, but worsens it in the case of the ESK model. Unfortunately, the x 2 /NDF ratio remains greater than one 
everywhere. It may be the consequence of our ignoring some important contributions, perhaps those with a scalar 
resonance considered in [4, 6]. The original models ESK, PB/HG, and AK do not contain any scalar resonances. We 
also have not included them as our main concern was the role of the cti resonance. It is also possible that the aipn 
Lagrangian should contain more terms than considered in (2.3). 

Using the D/S ratio as an additional data point is important. It can discriminate among various models. In our 
case it increases the separation of the ESK model from the others. It also strongly prefers larger values of the assumed 
a\ width. In calculations without the D/S ratio we were able to find a value of sin6> for each T ai that led to an 
acceptable fit to the e+e~ — > 27r + 27r~ cross section. However, the lowest value of r oi is excluded if we add the D/S 
ratio to the fitted data (see Table III). 

As to the partial decay width of p° — > 7r + 7r~7r + 7r~, the conclusion is not so categorical. Two models (PB/HG and 
AK) in their original forms provided results that were a little smaller, but did not contradict strongly the experimental 
value with its large errors. Only the original ESK model gave too large figure, which was in a clear disagreement 
with the experimental value. The inclusion of the a\iT intermediate states brought all values into the interval given 
by the experimental value and its errors summed linearly, see Table V. It must be said that the pure a\ model gives 
the decay widths that are close to the one-sigma upper limit or even beyond it. 

We originally hoped that our study would tell us the form of the aipn Lagrangian. But with respect to the aipn 
Lagrangian, no clear picture can be inferred from the low-energy results yet. The optimal values of the sine of the 
mixing angle, see Table IV, are from a broad interval and depend not only on the choice of the original model to which 
the ai diagrams are added, but also on the assumed value of the a\ width. The optimized values of sin 9 squeeze into 
interval (0.40,0.51). 

The quality of the fit over the whole energy range of the BaBar experiment [56], as measured by the x 2 /NDF ratio, 
seems to be better than that at low energies. But a more careful investigation in terms of the confidence level, which 
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TABLE VII: Survey of the mixing parameter sin 9 for various versions of two-component aipir Lagrangian (2.3) that appeared 
in literature. 



No. 


sin6< 


Reference 


1 





[38, 39] 


2 


0.2169 


[23, 24, 53] 


3 


0.5582 


[58] 


4 


0.6308 


[23, 24, 59] 


5 


1 


[37] 




0.40-0.51 
0.41-0.47 


our low-energy fits 
our all-energy fits 



takes into account \ 2 an d NDF separately, shows equal quality of those two fits. The values of sin# are compatible 
with those found at low energies, sin# S (0.41, 0.47) (lower boundary is obtained from the fits where T ai = 250 MeV 
was used). They occupy a narrower interval, what can be explained by lesser importance of sub-dominant diagrams 
with only p and 7r mesons in the intermediate states. 

It is interesting to compare our estimates of the mixing parameter sin#, defined in Eq. (2.3), with its values that 
have been used so far, see Table VII. The values in the first and fifth row simply reflect the fact that Lagrangians in 
Refs. [37-39] contained just one term. The remaining rows refer to various sets of four fundamental parameters (mo, 
g, a, and £) of the model in which the vector and axial-vector mesons were included as massive Yang-Mills fields of 
the SU(2)xSU(2) chiral symmetry [23]. The model built on previous works [28-30]. Our mixing parameter is related 
to the parameters r\\ and r\i of that model, which can be expressed in terms of the four fundamental parameters using 
Eqs. (2.9) and (2.10) in [23]. The formula is very simple 

• a ^ 2 

Vvi + V2 

In [23] , the fundamental parameters of the massive Yang-Mills model were determined using the experimental values of 
the masses and width of the p and a\ mesons. This procedure is not unique, there are two solutions. The corresponding 
mixing parameter is shown in the second and fourth row of Table VII. Row 3 corresponds to an ad hoc choice of 
fundamental parameters made in [58]. Last two rows show the range of our results obtained from various models 
(ESK, PB/HG, AK) and various assumed values of T 0l (250, 400, and 600 GeV). Unfortunately, there is no overlap 
with the rows above. The issue definitely requires more attention. The phcnomcnological models may be improved by 
including other intermediate states, perhaps those with scalar resonances. The parameters of the massive Yang-Mills 
model may be tuned by using richer experimental input (T(ai — > ^7) and D/S ratio as in [24, 32]) and more realistic 
formulas for relating theoretical parameters to experimental quantities, e.g., calculating the a\ width as a\ — > 3-zr 
instead of a\ — > pn. 

Our failure to obtain a more precise value of the mixing angle of the aipir Lagrangian suggests that it is necessary 
to make a simultaneous fit to data about several physical processes. The natural candidates are the e+e~ annihilation 
into various four-pion final states, the decay of the r lepton into neutrino and three or four pions, and the exclusive 
hadronic reactions of the type investigated, e.g., in [37]. 
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